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BULK  MODULUS  OF  ELASTICITY  OF  VARIOUS  ELASTOMERS:  THEORY  AND  EXPERIMENT 


INTRODUCTION 

The  mechanical  characterization  of  linear  homogeneous  isotropic  solids 
consists  of  the  identification  of  two  complex  elastic  moduli  as  functions  of 
frequency,  temperature,  and  pressure.  In  principle,  any  two  moduli  would  be 
suitable  since  any  pair  of  independent  moduli  may  be  converted  into  any  other 
pair  by  relations  well  known  from  elasticity  theory.  The  character  of  the 
experiments  designed  to  measure  a  particular  modulus  differs  considerably  from 
one  modulus  to  another,  and  the  ones  most  easily  measured  are  not  necessarily 
the  ones  needed  in  a  given  application  or  relevant  to  the  theory.  Problems 
also  arise  with  accuracy  in  converting  between  certain  pairs  of  elastic  moduli 
by  the  elastic  theory  relations  mentioned  above.  For  example,  measurements  of 
speed  of  sound  in  a  bar  is  a  favored  method  [1] .  The  corresponding  elastic 
modulus  so  obtained  is  Young’s  modulus  E.  This  value  may  be  combined  with  a 
torsional  or  other  determination  of  the  shear  modulus  G  [2,3] .  If  one  wishes 
to  obtain  either  the  bulk  modulus  K  or  the  plane-wave  (dilatational)  modulus  M 
from  E  and  G,  a  division  by  the  factor  (3G-E)  is  required.  However,  if 
Poisson  s  ratio  is  close  to  the  value  0.5,  as  it  is  for  most  elastomers,  the 
value  of  E  is  very  close  to  3G,  and  the  divisor  (3G-E)  becomes  very  small.  As 
a  consequence,  the  slightest  inaccuracies  in  either  G  or  E  are  magnified 
greatly,  and  the  resulting  values  of  K  or  M  have  such  large  uncertainties  as 
to  mask  any  finer  detail  that  might  be  of  significance.  For  this  and  other 
reasons,  it  is  desirable  to  have  a  direct  measurement  of  the  bulk  modulus  K 
or  the  plane-wave  modulus  M  [4] .  ’ 

.  ,  ' Dne,°f  thes?  othfr  reasons  is  that,  from  a  theoretical  standpoint,  the 
bulk  modulus  arises  from  quite  different  mechanisms  than  the  other  elastic 
mo  u  i  i  particular ,  it  is  mostly  a  measure  of  the  short-range  inter-  and 
intramolecular  potentials  and  is  little  affected  by  longer  range  chain  effects 
that  dominate  m  shear  and  tension.  Thus  the  theoretical  information  on 

d^r^lar  S?ale  “echaaicsL furnished  by  the  bulk  modulus  is  to  an  important 
degree  complementary  to  that  obtained  from  the  other  moduli.  The  bulk  moduli 
of  elastomers  above  their  glass-transition  temperatures  are  quite  similar  to 
those  of  organic  liquids  both  in  magnitude  and  in  loss  characteristics.  This 

dL^"1“te“ii.W1khiihe  id?a  that  short-ranSe  interactions  dominate  in 
determining  the  bulk  moduli  of  both  liquids  and  elastomers. 

UndeJwatereSound°R  ?f  &  ^  measure  dynamic  bulk  moduli  at  the 

USRDW-!  ?  i  Re^erence  Detachment  of  the  Naval  Research  Laboratory  (NRL- 

caSbilitv  to°  UP°?,th!  d6?ire  to  have  an  in-house  experimental7  ( 

toPbe  described^n^V11  ***  d®velopmenfc  of  new  elastomers.  Since  the  method 
frequencies  easi1^  be  extended  to  quite  low  audio 

g  s  moduli  can  be  compared  at  the  same  frequency. 
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The  bulk  modulus  is  measured  at  NRL-USRD  in  an  acoustic  coupler  similar  to 
those  used  for  reciprocity  calibration  of  hydrophones.  The  method  is  very 
close  to  that  of  McKinney,  Edelman,  and  Marvin  [5] .  The  method  will  be 
described,  the  basic  equations  given,  calibration  and  measurement  techniques 
discussed,  and  preliminary  results  given  for  a  number  of  different  elastomers. 
In  addition,  results  of  a  detailed  analysis  will  be  given  for  several 
elastomers  to  show  the  capability  of  the  method  to  provide  comparatively 
complete  information  about  the  bulk  modulus  of  a  material  as  a  function  of 
both  temperature  and  frequency. 

THEORY 

Basic  Free-Volume  Model 

Because  of  the  close  resemblance  between  the  bulk  moduli  of  elastomers  and 
typical  organic  liquids,  it  appears  likely  that  the  free-volume  theory  of 
liquids  can  provide  a  suitable  framework  within  which  to  understand  the  bulk 
moduli  of  the  elastomers  of  interest  in  the  present  work. 

In  its  basic  form,  the  free-volume  concept  envisions  a  liquid  above  its 
glass-transition  temperature  as  a  somewhat  spongy  material  consisting  of  many 
holes  of  molecular  size  imbedded  in  a  highly  compliant,  loose  structure. 

These  holes  vary  widely  in  size  but  are  distributed  randomly  throughout  the 
material  at  any  given  instant.  At  or  below  the  glass  temperature  T  ,  the 

material  is  in  the  form  of  an  amorphous,  glassy  solid.  Such  solids  are  quasi- 
crystalline  in  structure  with  the  molecules  almost,  but  not  quite,  arranged  in 
the  regular  orderly  manner  of  a  crystalline  solid.  About  each  molecule  in  the 
glassy  state  there  is  a  local  volume  over  and  above  that  actually  occupied  by 
the  molecule  itself.  This  local  volume  v^  is  the  "free  volume."  It  forms 

"cages"  on  a  molecular  scale  within  which  each  molecule  is  constrained  to  move 
at  low  temperatures.  There  is  little  or  no  molecular  transport  possible  below 
the  glass-transition  temperature. 

At  Tg  a  phase  transition  takes  place,  a  sort  of  "melting"  of  the  glassy 

state;  this  phase  change  is  accompanied  by  a  marked  increase  in  the  thermal 
expansion  coefficient,  a^,,  of  the  material,  by  more  than  two  orders  of 

magnitude,  as  a  rule.  As  the  temperature  rises  further,  the  excess  volume 
that  accompanies  the  thermal  expansion  is  distributed  throughout  the  sample, 
but  not  uniformly.  The  excess  free  volume  is,  instead,  randomly  distributed 
among  the  N  molecules  of  the  sample,  and  some  molecules  therefore  have 
relatively  large  holes  around  them,  while  others  have  only  small  holes  as 
their  share  of  the  excess  free  volume.  This  process  of  redistribution  of  free 
volume  is  a  continually  changing  one,  but  the  time  scale  for  such  changes  is 

estimated  to  be  of  the  order  of  lO-12  second  in  typical  liquids  [6],  long 
enough  on  the  time  scale  of  molecular  motions  for  one  to  adopt  a  quasi— static 
picture  of  the  structure  within  which,  for  the  duration  of  its  existence,  one 
can  calculate  the  distribution  of  free-volume  hole  sizes.  Cohen  and  Turnbull 
[7]  have  given  a  derivation  of  the  distribution  function  for  hole  sizes  by 
considering  the  excess  total  free-volume  to  consist  of  packets,  one  for 

each  of  the  molecules  in  the  sample  with  no  restriction  that  the  packets  all 
be  of  the  same  size.  A  simpler  derivation,  restricted  to  packets  of  the  same 
size,  has  been  given  by  Bueche  [8]  (see  also  Appendix  A).  If  all  packets  are 
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of  the  same  size,  it  follows  that  the  chance  for  a  particular  molecule,  or 
molecular  segment  in  the  case  of  a  polymer  chain,  to  have  q=vl./v-  packets 
associated  with  it  is 

.  P(vjs)  =  (Vf^Tvp1/2  exp(-avj,/vf)  ,  (1) 

where  a=ln  q-1.  Turnbull  and  Cohen  [9]  and  more  recently  Cohen  and  Grefct  [6] 
defined  a  critical  value  v^  for  the  excess  free  volume  such  that  holes  larger 
than  Vj  can  contain  molecules  or  molecular  segments  that  are  able  to  move 
about  freely  within  the  hole  without  encountering  any  energy  barriers,  much  as 
a  gas  molecule  might  move  within  v^.  We  adopt  this  two-phase,  gas- 
phase  model  for  the  rubbery  state  of  an  elastomer  above  the  glass  temperature. 

The  elastomer  is  considered  to  contain  holes  of  varying  sizes,  some  of 

which  are  larger  than  v^  and  may  contain  free-moving  segments  of  molecules 

that  behave  like  the  molecules  of  a  gas.  The  total  pressure  of  such  a  gas 
phase  in  a  sample  at  temperature  T  may  be  written  as  for  an  ideal  gas: 


P  =  N*kT/vJ  (2) 

■where  N  is  the  total  number  of  molecules  (or  segments)  that  are  in  the  gas 

phase  in  holes  larger  than  vt.  The  total  volume  of  these  large  holes  is 
*  * 

V*. 

We  note  that  the  isothermal  bulk  modulus  ICp  for  an  ideal  gas  is  equal  to 

its  pressure,  and  the  adiabatic  modulus  is  5/3  the  isothermal  modulus. 

This  relates  Eq.  (2)  to  the  bulk  modulus;  our  model  is  based  on  this  equality. 
The  model  postulates  that,  as  far  as  the  bulk  modulus  is  concerned,  the 
elastomer  behaves  like  an  ideal  gas  contained  in  holes  within  a  highly 
compliant  matrix  so  the  external  pressure  is  transmitted  freely  throughout  the 
solid,  as  it  is  in  the  case  of  a  liquid.  Corrections  for  the  residual 
rigidity  of  the  matrix  can  be  made  if  necessary.  They  may  be  incorporated  by 
considering  the  compliance  of  the  elastomer  to  be  the  sum  of  the  compliances 
of  the  gas  and  the  matrix;  i.e.,  the  compliances  are  "in  series." 

can  be  calculated  from  the  hole  size  function,  Eq.  (1) .  The  result  is 


Vf  =  A  exp  (-aNv*/Vf)  ,  (3) 

where  A  and  a  are  constants,  N  is  the  number  of  molecules  in  the  sample,  and 
“f  is  the  total  sample  volume.  Note  that  V^Nv. . 
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In  order  for  a  molecule  to  leave  the  "wall"  of  a  hole  and  enter  the  gas 
phase,  it  must  free  itself  from  the  forces  that  bind  it  to  the  wall,  namely 
the  attractive  forces  of  its  neighboring  wall  molecules.  This  evaporative 
process  requires  energy,  and  it  follows  that  N  must  be  proportional  to  a 
Boltzmann  factor  of  form  exp(-e/kT)  where  e  is  the  energy  of  vaporization. 
This  energy  is  not  a  constant  but  depends  upon  the  hole  size  for  the  same 
reason  that  the  vapor  pressure  in  a  small  bubble  depends  upon  the  size  of  the 
bubble.  As  the  free  volume  vanishes  at  the  glass-transition  temperature,  e 
must  also  go  to  zero.  In  large  holes  that  can  exist  at  high  temperatures,  e 
approaches  a  value  comparable  to  the  energy  of  vaporization  of  the  material. 

Putting  these  factors  together  into  Eq.  (2)  gives,  finally, 


K  =  AT  exp(C/Vf-e/kT) . 


(4) 


where  C=aNv£. 

The  total  free-volume  in  Eq.  (4)  is  proportional  to  the  total  volume  V 

of  the  sample  and,  therefore,  it  varies  linearly  with  temperature  and  external 
pressure  as 


vf  - 


fg 


0CTV  (T-T  )  -  P  P 
T  gv  g'  rs 


(5) 


where  a^,  is  an  appropriate  volume  thermal  expansion  coefficient,  V  is  the 

total  volume  at  T  ,  and  /?  is  the  compressibility  of  the  solid.  Since  6  =K-1 
S  s  '  s  y 

inclusion  of  the  p  P  term  makes  Eq.  (5)  for  K  inherently  nonlinear.  However, 
at  modest  pressures,  the  /JgP  term  is  small  and  an  average  value  of  Ps  may  be 
used  without  appreciable  error.  This  linearizes  Eq.  (5).  In  what  follows,  we 
take  P  0  and  omit  the  ps P  term,  treating  only  the  temperature  and  frequency 
dependence  of  K.  Extension  to  higher  pressures  is  straightforward. 


V-  at 
fg 


In  Eq.  (5),  V,  is  defined  in  such  a  way  that  it  has  the  value  ,  „  x 

,  1  fg  g 

so  the  correct  expansion  coefficient  to  use  in  Eq.  (5)  is  the  difference 

between  the  expansion  coefficients  for  the  liquid  phase  and  the  extrapolation 
of  the  coefficient  of  the  glassy  phase.  Figure  1  below  shows  more  clearly  the 
meaning  of  this  statement. 
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Fig.  1  -  Thermal  expansion  coefficients  of  a  typical  elastomer  in  the 
vicinity  of  the  glass-transition  temperature. 


The  value  of  ccT  to  be  used  in  Eq.  (5)  is  a^g-otg,  and  the  temperature 
variation  of  is  given  by  Eq.  (5)  with  this  value  of 

Substituting  Eq.  (5)  into  Eq.  (4)  and  taking  the  logarithm  of  the 
resulting  expression  gives 


£n(K/T)  =  A'  -  B/T  + 


Vf  +  “TVg  (T'TJ 


ln[VaT  YT-yj 


(6) 


where  A'=Ngk/2,  C=aNv*,  Vfg,  and  aT  are  all  constants  for  a  given  material. 

Ns  is  the  number  of  "surface"  molecules  or  molecular  segments  lining  all  holes 
with  vf>vf.  The  first  term  on  the  right  is  negligible  compared  with  the  third 
because  C/Nvf  *vj/vf  ^40,  and  A1 /N^Nsk/N«0.025.  Also,  the  logarithmic  term 
is  small  compared  with  the  C/Vf  term  so  Eq.  (6)  may  be  written  more  simply 


In  (K/T) 


+  otmV  (T-T  ) 
fg  T  gv  gJ 


(7) 


a  g0°V aPPro^imation  not  too  far  above  the  glass-transition 
linear  form'  ”  convenient  for  curve  fitting  to  write  this  in  the 


[in  (K/T)] 


=  Vfg  +  aTVg 


(T-Tg)  =  d  + 


bT 


(8) 
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WLF  Frequency  Shift  Property  of  the  Model 

It  is  a  simple  matter  to  show  that  Eq.  (7)  leads  to  the  Williams,  Landel, 
Ferry  [11]  equation  for  the  shift  in  frequency  required  to  superimpose  two 
values  of  K1  or  K"  measured  at  different  temperatures.  This  is  the  familiar 
WLF  frequency  shift  method,  though  it  will  not  appear  here  in  its  conventional 
form  since  we  have  not  yet  introduced  the  frequency  dependence  of  K. 

Consider  two  temperatures,  T^  and  Tg .  Write  the  free  volume  at  T2  as 

vf2  =  Vf,  *  “TVl  <W  (8) 


where  V  is  the  total  volume  at  T  .  Then,  neglecting  second  order  terms  in 
a,p,  Eq.  (7)  becomes 


W./W  -  .  (10) 

where 

c1=C/Vfl  and  c2  =  aT_1  V£l/Vr  (11) 

We  note  that  the  right-hand  side  of  Eq.  (10)  has  the  form  of  the  WLF 
equation  for  ln(a^,),  the  logarithmic  frequency  shift  factor  to  make  data  taken 

at  T  correspond  with  those  taken  at  T„  [12] .  The  constants  c,  and  c  are 

identified  with  the  material  properties  of  the  sample  through  Eq.  (11) * 

However,  c^  contains  which  depends  upon  the  particular  choice  of  reference 

temperature  T^  .  Bueche  [10]  has  suggested  that  the  product  ^^=2080  for  a 

large  number  of  materials,  and  Ferry  [13]  has  also  shown  that  the  product  c  c 
is  a  constant  independent  of  the  reference  temperature.  Our  data  do  not  *  2 
support  this  universal  value  hypothesis;  our  values  of  c  and  c  and  their 
product  vary  appreciably  from  the  "universal"  values.  1 


High-Temperature  Behavior  and  the  Arrhenius  Factor 


Before  introducing  dynamic  factors  we  note  that,  at  low  temperatures  below 
about  100  to  150° C  above  the  glass  temperature  T^,  the  e/kT  term  is  small 

compared  with  the  other  term,  C/V^,  so  the  latter  dominates  at  lower 


temperatures.  However,  at  temperatures  above  about  Ty+100°C,  e/kT  becomes 

large  enough  to  influence  the  modulus  and,  as  found  experimentally,  this  leads 
at  higher  temperatures  to  an  "Arrhenius"  type  of  behavior  of  K  with 
temperature  due  to  the  exp  (-e/kT)  term.  We  shall  disregard  the  Arrhenius  term 
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in  what  follows.  Nevertheless,  the  present  model  introduces  this  high- 
temperature  behavior  in  a  straightforward  way  and  shows  how  the  transition 
takes  place  between  low-  and  high-temperature  behavior. 

Dynamical  Considerations  and  Frequency  Dependence 

The  frequency  dependence  of  K  may  be  introduced  by  way  of  relaxation 
theory  which  is  fundamentally  related  to  causality  and  is  a  necessary  property 
of  all  macroscopic  linear  systems.  Briefly,  a  linear  system  with  a  single 
relaxation  time  r  responds  to  an  input  F(t)  to  give  an  output  of  the  form 


H(t) 


ft 

F(t')  exp [(t1 -t) /t]  dt 1  , 


o 


(12) 


where  F(t' )=0  for  all  times  t’<0.  This  simply  means  that  the  output  at  any 
time  is  the  sum  of  all  inputs  at  earlier  times,  each  declining  exponentially 
in  amplitude  with  time.  Applied  to  a  sinusoidal  input  of  circular  frequency, 
w,  the  result  expressed  in  terms  of  the  complex  bulk  modulus,  K*,  takes  the 
well-known  form  [14] 


K*  =  K' 

+  iK" 

(13) 

K'  =  K' 
o 

♦  4K— - 

1  +  (wr)2 

(14) 

K"  =  AK 

WT 

(15) 

1  +  (wr)2 

where  is  the  limiting  value  of  K'  as  frequency  goes  to  zero,  and  AltK^-K' 

where  is  the  limiting  value  at  high  frequency. 


In  the  equations  above,  the  complex  bulk  modulus  K* 
an  in-phase  "storage  modulus,"  K' ,  and  an  out-of-phase 


is  the  combination  of 
"loss  modulus,"  K" . 


Define  a  logarithmic  frequency  scale  by 


X  -  ln(f/fm)  where  =  1/2 jtt  and  f  =  u/2 T  . 


Then  Eqs.  (14)  and  (15)  can  be  written  in  the  form 


K'  =  k;  +  f  (1+tanh  x) 


(16) 


(17) 
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K"  =  ~  sech  x  . 


(18) 


In  their  study  of  dielectric  losses  in  polymers,  Fuoss  and  Kirkwood  [15] 
showed  that  the  imaginary  part  of  the  dielectric  constant  that  has  the  same 
form  as  Eq.  (18)  could  be  generalized  to  the  case  of  a  continuous  distribution 
of  relaxation  times  by  simply  introducing  a  "spread  factor,"  a,  that 
multiplies  x  in  Eq.  (18).  (cc  .  is  positive,  real,  and  <  1.)  They  also  showed 
that  the  spread  factor  determines  a  relaxation  time  spectrum  of  the  polymer. 

We  adopt  this  generalization  for  the  bulk  modulus  and  write 


K'  =  Ko  +  Y  [1+tanh  (ax)]  (19) 


K"  =  2~  sech  (ax)  •  (20) 

Appendix  B  gives  an  abbreviated  form  of  the  Fuoss  and  Kirkwood 
generalization  and  their  expression  for  the  relaxation  time  spectrum,  G(r), 
which  depends  upon  a . 

If  the  principle  underlying  the  WLF  shift  method  is  to  hold,  it  is 
necessary  that  curves  of  K'  and  K"  vs  frequency  measured  at  different 
temperatures  be  similar  in  shape  so  they  can  be  superimposed  by  suitable 
shifting. along  the  frequency  scale,  the  amount  of  such  shift  being  determined 
by  the  difference  in  temperature  of  a  given  curve  relative  to  a  chosen 
reference  temperature.  It  is  therefore  required  that  K' (T  )=K"(T  )  under  a 
frequency  shift  from  x(T  )  to  x(T  )  where  1  ° 


*(TC)  =  *(T.)  -  fc.-x.  )  . 


(21) 


The . subscript  "o"  identifies  a  reference  temperature  which  may  be  chosen 
arbitrarily  and  the  subscript  i  identifies  the  inflection  point  of  the  K1 (x) 
curve  or  the  maximum  of  the  K" (x)  curve. 

Figure  2  shows  the  rationale  for  Eq.  (21). 
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It  is  clear  from  Fig.  2  that  a  log  frequency  shift  of  the  amount  (x.  -x.  ) 
is  needed  to  superimpose  the  curve  upon  the  Tq  one.  Thus,  any  point ""x  on° 
the  Tt  becomes  a  corresponding  point  on  the  Tq  curve  under  the  shift 

(xil  ~xio) »  and  E(l-  (21)  expresses  this.  Since  x=ln(f/fm) ,  Eq.  (21)  may  be 
written 

X(Ti  =  x(V  -  (xil-xio>'  or  lntf(T0)/f(Tl)]  =  <22) 

Thus, 


It  follows  that  (x-xm)  is  invariant  with  respect  to  changes  in 


(23) 


1U  *  - - O  ~  ~ 

temperature,  and. this  is  the  proper  variable  to  use  in  the  argument  of  the 
written110  functlons  ln  K'  and  K" *  Therefore,  Eqs.  (19)  and  (20)  should  be 


K'  “  Ko  +  T  {1+hanh  [a(x-xm)]> 
K"  =  |^  sech  [a(x-xm)J 

with 


(24) 

(25) 
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x  (T)  -  x  (T  )  = 

mv  7  mv  o7 


(26) 


where  c.  is  the  value  of 


c.  at  reference  temperature  T  . 


J  o 


The  frequency  dependence  of  the  model  is  contained  in  Eqs.  (24)  and  (25) 
which,  as  we  have  shown,  also  exhibit  the  WLF  shift  property.  At  present,  we 
are  uncertain  how  these  equations  should  be  combined  with  the  zero  frequency 
expression  Eq.  (4)  to  produce  a  complete  description  of  both  frequency  and 
temperature  dependence  in  a  single  expression.  A  simple  approach  would  be  to 
argue  that  since  the  loss  modulus  KB  is  proportional  to  the  viscosity,  the  AK 
factor  that  multiplies  sech(x-xm)  should  have  the  form  of  Eq.  (4). 

The  losses  depend  upon  the  viscous  forces  acting  in  the  medium,  whether 
solid  or  liquid.  In  fact,  the  lossy  part  of  the  Young’s  and  shear  moduli  are 
proportional  to  the  coefficient  of  shear  viscosity,  fj  [16,17].  It  is 
reasonable  to  assume  the  existence  of  a  compressional  viscosity  which 
determines  the  loss  bulk  modulus  K"  in  the  same  way,  although  the  loss 
mechanisms  will  not  be  the  same.  In  the  shear  case,  the  viscosities  of  many 
organic  liquids  obey  a  relation  of  the  form 


V'  (T)  =  T)>  (To)|-  exp  [C/Vf  (T)-C/Vf  (To)]  . 


(27) 


o 


This  was  determined  empirically  by  Doolittle  [18]  and  is  commonly  referred  to 
as  the  Doolittle  equation.  In  Eq.  (27)  it  has  been  cast  in  the  form  suggested 
by  the  free-volume  theory  of  liquids  [19]  ,  and  fj  is  proportional  to  one  of  the 
factors  in  Eq.  (4) .  Thus  we  already  have  a  viscosity  term  in  the  expression 
for  K  at  very  low  frequency.  Since  Kn  must  be  proportional  to  the  viscosity, 
it  follows  that  K"  in  Eq.  (25)  must  be  modified  to  incorporate  the  viscosity 
term.  This  modification  is  guided  by  the  condition  that,  as  frequency  goes  to 
zero,  the  bulk  modulus  must  reduce  to  the  form  of  Eq.  (4).  The  final  result 
is  that  AK  is  proportional  to  the  right-hand  side  of  Eq.  (4)  which  makes  AK  a 
function  of  T  instead  of  a  constant,  as  it  is  m  Eqs.  (24)  and  (25) .  This,  in 
turn,  destroys  the  frequency  shift  property  because  the  frequency  shift  in 
this  case  does  more  than  simply  move  the  K'  and  K”  curves  along  the  frequency 
axis  unaltered  in  shape;  a  frequency  shift  is  also  accompanied  by  a  change  of 
shape  through  the  temperature  dependence  of  AK. 

For  the  present,  we  will  treat  the  frequency  dependence  as  given  by  Eqs. 
(24)  and  (25)  in  an  ad  hoc  manner,  considering  Eq.  (4)  to  be  the  fundamental 
result  of  the  free-volume  approach  via  our  model.  We  therefore  take  AK  to  be 
a  constant  for  a  given  sample  and  use  the  WLF  method  to  construct  master 
curves  of  K1  and  Kn  as  functions  of  frequency  at  a  reference  temperature  T 

° 

This  process  determines  the  parameters  c{  and  ^  of  the  WLF  shift  factor  3^ 
where,  by  definition,  log  a^  =  -(x-xj,  or 
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ln(aT)  =  -(x-xm)  = 


~C1  (T-Tq) 

c2  +  (T-T0)  ’ 


(28) 


In  addition  to  these  parameters,  the  spread  factor  a, 

a/6  ^ounc*  by  fitting  the  curve  obtained  by  the 
of  Eqs .  (24)  and  (25).  This  completes  the  parametric 
master  curves. 


and  the  parameters  K' 

o 

WLF  method  to  the  forms 
description  of  the 


be  difficult6  ^  aU>  ^  fittinS  «  the  general  case  proves  to 

51*2  •  Fortunately,  in  many  cases,  some  of  these  parameters  can  be 

mine  y  methods  other  than  a  brute  force  computer  approach.  For 
example,  if  the  data  at  the  reference  temperature  Tq  include  the  inflection 

point  of  K' ,  the  value  of  fffl  can  be  found  by  inspection.  Differentiation  of 
Eq.  (24)  with  respect  to  x  gives 


dK'/dx  =  a~  sech2  [a(x-xm)] 


(29) 


sl°Pes  to  *  hyperbolic  secant-squared 

WTF  f  f  of  Eq-  (29)  also  determines  the  parameters  a,  AK,  and  the 

WLF  parameters  and  c2 .  (See  Appendix  C)  ’  6 

Ferry  [19]  has  shown  that  the  product  c2  should  be  a  constant 

independent  of  temperature.  This  means  that  only  one  of  these  c’s  is  an 

neededntoDfir^etr'  ^  ^  °nly  three  P^^ers  would  be 

“®®,Jl  t0,  ,  ^he.?l?P®s  to  Eq-  (29),  and  the  fitting  process  would  be 

° 51  r:b  1 -ddif?ur  df a  gera\iy  d° 

1  2  constant  for  different  samples,  however,  and  this  is  puzzling  since 

it  runs  counter  to  predictions  of  the  WLF  theory.  Ferry  himself  [19]  gives 
experimental  data  on  shear  moduli  together  with  values  of  c  and  c  which  do 
not  fully  satisfy  the  constant  product  requirement  of  theory. 

tha^th^fractiona^fiee  iolLr^the^lasi"1 1  .  othfrs’  have  suggested 

coefficient  is  approximately  4.8X10H  CO"1  for  most  liquids  in  the 
temperature  range  up  to  about  100*0  above  T  With  these  values  of  free 

volume  and  expansion  coefficient  thp  WTF  *  +.  , 

„  .  cieM»  tbe  WL*  ^  factor  may  be  written  in 

universal"  form 


11 


BURNS,  DUBBELDAY,  and  TING 


-40  (T-T  ) 

ln(aT)  =  52  +  (T-T  )  •  (30) 

g' 

This  is  the  form  suggested  by  Bueche.  Parenthetically,  it  implies  that  ~ 
4°Vfg’  that  is’  a  hole  about  40  times  larger  than  the  average  size  of  a  hole 
s-t  Tg  is  needed  to  accommodate  a  molecular  segment.  Thus,  many  n average** 

holes  must  aggregate  in  order  to  form  a  hole  of  sufficient  size  into  which  a 
segment  can  jump  in  the  liquid  or  rubbery  state.  Bueche’ s  argument  shows  that 
C1  anc*  c2  are  rebate(i  to  physical  constants  of  a  given  liquid  and  by 

inference  to  the  corresponding  constants  of  a  given  elastomer.  Unless  these 
constants  vary  more  widely  in  elastomers  than  in  liquids,  which  is  certainly 
possible,  and,  unless  they  are  much  less  interrelated  in  elastomers  than  in 
liquids,  it  is  surprising  that  we  do  not  find  their  product  to  be  more  nearly 
constant  in  the  materials  we  have  studied.  As  mentioned  above,  the  data  Ferry 
gives  [19]  on  a  variety  of  different  elastomers  also  fail  to  support  the 
constancy  of  the  product  c^ c^ .  No  reason  for  this  has  been  suggested. 


Relaxation  Spectrum 

From  a  theoretical  point  of  view,  one  of  the  more  interesting  kinds  of 
information  one  would  like  to  get  from  elastic  modulus  measurements  on 
elastomers  is  the  relaxation-time  spectrum  associated  with  the  loss  mechanisms 
in  the  materials. 

Fuoss  and  Kirkwood  [15]  have  shown  that  the  spread  factor  determines  a 
relaxation  spectrum  given  by 


G(r) 


OCT 


XT 


cos(gx/2)  coshCgsl 


cos2 (ax/2)  +  sinh2 (gs) 


(31) 


where  s  is  defined  as  -ln(r/rffl) .  (See  Appendix  B  for  details.) 

It  should  be  pointed  out  that  only  two  parameters  are  required  to 
completely  determine  the  relaxation  spectrum:  the  spread  factor  and  the  time- 

scale  factor  Tm  or  its  equivalent  fffl (==7"m  *),  the  frequency  at  the  inflection 

point,  of  K  or  at  the  maximum  of  K  .  The  spectrum  G (r)  is  a  1/r  spectrum 
modified  by  the  factor  in  brackets  in  Eq.  (31) .  When  plotted  against  the 
logarithm  of  T,  Eq.  (31)  has  the  form  shown  in  Fig.  3  below.  The  result 
clearly  depends  on  the  assigned  value  for  the  spread  factor  a. 


12 


NRL  MEMORANDUM  REPORT  NO.  5991 


In  general,  this  type  of  spectrum  is  typical  of  liquid  behavior.  The 
experimental  value  of  fm  falls  in  the  range  0.1  to  1  kHz  for  several  elastomer 

samples  measured  at  NRL-USRD  so  rm  is  about  1  millisecond,  and  the  spectrum 

typically  extends  to  relaxation  times  well  below  0.1  nsec  with  short 

rh^a£10n  £lme?  °f^he  order  of  1  /isec  or  less  predominating.  The  1/r 
samn^Cter  of  relaxation  spectra  exhibited  experimentally  by  our  elastomer 
mples  are  typical  of  liquids  and  are  consistent  with  free-volume  theories 
for  normal  liquids  (Cohen  and  Crest  [6]).  By  contrast,  various  chain  Seories 

[231  °fThendiffer^  ^  m°duli  ^edict  more  or  less  r"1 /2  relaxation  spectra 
L  J  •  difference  in  r-dependence  between  shear  and  bulk  (liquid-)  theories 

is  pronounced  and  is  an  indicator  of  the  difference  in  molecular  mechanisms 
that  operate  in  bulk  as  compared  with  shear  or  tension.  In  lipids  short 

in^hear  itais10nf  ar?  do”inant  in  determining  compressive  properties, ^hile 
primarily  longer  range  chain  entanglements  and  interactions 
that  gave  r1Se  to  energy  transfers  that  accent  for  relaxation  p^oceLe 

Fuosetand°Kirk»ood°if  adprediction‘of  their  th  ““  rex“ation  sP«tra  given  by 

experimentally.  Therefore  it  ms  anc*  tas  not  been  verified 

sSch  Sp,^ediitL^KnHoSv:“  - in  ”lyins  - 

expect  Wi^a  rerpoSiororK^^dV''6"^3  in  *  Mterial  ««!■» 
parameters:  a,  AK  K*  and  f  Jj,*  K  c?rve?’  e*ch  with  a  set  of 

’  ’  o’  and  fm‘  There  would  also  be  a  superposition  of 
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relaxation  spectra  G(r)  centered  on  different  r  and  each  with  a  different  a 

m 

Such  superpositions  can  make  the  experimental  K1  and  K"  curves  appear  quite 
different  from  the  simple  form  we  have  dealt  with  so  far,  particularly  if  the 
spectra  of  the  different  relaxation  processes  are  not  widely  separated  in 
frequency.  We  note  this  generalization  but  will  not  deal  with  it  further  in 
this  report.  We  have  not  seen  clear  evidence  for  more  complex  K1  curves  in 
the  limited  experimental  data  obtained  to  date. 


EXPERIMENTAL 

There  are  two  principal  ways  in  which  the  dynamic  bulk  modulus  of  a  solid 
material  can  be  found  experimentally.  One  is  the  impedance  tube  method  in 
which  the  speed  of  transmission  of  a  plane  wave  through  a  plane  sample  is 
measured.  From  this,  the  plane-wave  modulus  M  is  found,  and  the  bulk  modulus 
K  can  be  calculated  from  the  fundamental  relation  K=M-4G/3,  provided  the  shear 
modulus  G  is  also  known.  As  mentioned  earlier,  in  the  materials  of  interest 
in  the  present  work  G  is  at  least  two  orders  of  magnitude  smaller  than  M,  so  K 
equals  M  to  within  1%  or  better.  There  are  technical  difficulties  in  making 
accurate  measurements  by  the  impedance  tube  method,  mostly  connected  with  edge 
effects  when  the  sound  wavelength  is  not  negligibly  small  compared  with  the 
dimensions  of  the  sample.  However,  useful  results  can  be  obtained  with  this 
method  with  due  care.  It  is  especially  useful  at  frequencies  in  the  kHz  range 
and  complements  the  acoustic  coupler  method  which  is  best  at  lower 
frequencies . 

The  second  method  is  the  acoustic  coupler  method  developed  in  1955  at  the 
National  Bureau  of  Standards  by  McKinney,  Edelman,  and  Marvin  [5].  The  method 
employs  a  heavy-wall,  steel  pressure  vessel  in  which  there  is  a  cavity 
containing  two  piezoelectric  transducers,  in  addition  to  the  sample  and  a 
filling  fluid  whose  compressibility  is  accurately  known  as  a  function  of  both 
temperature  and  pressure.  One  of  the  transducers,  the  "driver,"  acts  as  a 
volume  expander  and,  as  its  volume  changes  in  response  to  an  applied  voltage, 
a  pressure  change  is  produced  in  the  chamber.  The  pressure  acts  on  the 
receiver  transducer  to  produce  an  output  voltage  that  is  a  measure  of  the 
resulting  pressure  change  in  the  chamber.  Therefore,  the  ratio  of  output  to 
input  transducer  voltages  carries  information  about  the  compliance  of  the 
coupler  contents  and  particularly  about  the  compliance  of  the  sample. 

The  coupler  used  at  NRL-USRD  for  the  calibration  of  hydrophones  has  been 
described  by  Bobber  [24] .  For  the  present  application  the  same  kind  of  heavy 
steel  chamber  is  used  except  that  only  two  of  the  customary  three  ports  are 
needed,  one  for  a  driver  transducer  and  one  for  a  receiver.  Reciprocity  is 
not  required  of  the  transducers  in  the  present  case.  To  make  this  point  clear 
and  to  show  the  approximations  made  in  the  basic  equations,  the  main  line  of 
the  derivation  is  given  here.  It  follows  closely  the  analysis  of  McKinney, 
hdeiman,  and  Marvin. 

The  sum  of  the  changes  of  volume  AVi  of  all  the  items  inside  the  coupler 
equals  the  change  of  volume  of  the  cavity,  AY  ,  or 
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AY  =  I  AV.  . 

c  .  1 

1 


(32) 


It  is  essential  to  assume  that  the  pressure  changes  AP  throughout  the  cavity 
are  uniform;  i.e.,  AP  must  be  the  same  everywhere  in  the  cavity  during  the 
pressure  cycle  produced  by  the  driver  which  acts  as  a  volume  expander.  This 
sets  an  upper  limit  to  the  frequency  that  may  be  employed  in  a  given  chamber. 
Dividing  Eq.  (32)  by  AP, 


AV  r  AY. 

c  _  y  i 

AP  .  AP  (33) 

i 

For  the  passive  material;  i.e. ,  everything  in  the  chamber  including  the 
chamber  itself  but  excluding  the  transducers,  one  may  set 


1  AV 
V  AP 


1  dV 
V  dP 


=  P  =  K 


-1 


(34) 


where  K  is  the  bulk  modulus  and  /J  is  the  compressibility.  For  active 
materials,  specifically  the  piezoelectric  transducers  in  the  chamber,  one  may 
write  in  the  linear  region,  J 


E  =  aP  +  bu 


and 


(35) 


I  =  cP  +  du 


where  E  is  the  voltage  across  the  transducer,  I  is  the  current  through  the 
transducer,  u  is  the  volume  velocity,  P  is  the  pressure  in  the  cavity,  and  a 
;n,d  d  are  properties  of  the  transducer  material.  These  relations  are  ’ 

2  resicti)elv  l  “  V*  C°"pler-  indic!“ed  ibices  1  and 

*  n  ls.not  necessary  to  assume  reciprocity  of  either 

iranSrJUCer*  ^  algeknaic  manipulation,  the  electrical  parameters  of  the  two 
transducers  may  be  related  to  the  respective  material  parameters  of  the  driver 

Kr  F°-  lnSt?"Ce!  *»«■>  of  driver  to  receiver  voltages  is  £ve£ 

for  the  open  receiver  circuit  case  by  6  given 


E1 

R  =  r  = 


bid2 


C2  ai 

iwC  -  /  - 

 d2  blj 


b  c  -ad 

2  2  2  2 


(36) 
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where  C  is  the  total  compliance  of  all  the  passive  materials  in  the  coupler 
minus  the  chamber  compliance. 

To  calibrate  the  chamber,  one  measures  the  above  voltage  ratio  Rq  when  the 

chamber  is  filled  with  a  fluid  of  known  compressibility  /?q,  and  again  when  the 

coupler  contains  a  metal  of  known  volume  and  known  (small)  compressibility 

In  this  case,  the  voltage  ratio  is  R^.  Measuring  the  voltage  ratio  Rg 

with  an  elastomer  sample  of  volume  Vg  in  the  coupler,  the  compressibility  p 

of  the  sample  may  be  calculated  from  the  following  expression  which  is  easily 
derived  from  Eq.  (36) .  It  gives  the  sample  compressibility  in  terms  of 
measured  voltage  ratios  for  the  coupler  containing,  variously:  brass,  the 
unknown  sample,  and  only  oil.  The  relation  is 


(37) 


All  of  the  R  ratios  in  this  expression  are  complex  numbers,  as  Eq.  (36) 
suggests.  In  general,  p '  will  be  a  complex  number  whose  real  part  gives  the 

storage  compressibility  of  the  sample  and  whose  imaginary  part  gives  the  loss 

compliance.  Both  p q  and  p ^  are,  for  all  practical  purposes,  real  because  the 

fluids  used  in  these  measurements  and  the  metal  calibration  materials  (brass 
in  the  present  case)  have  very  low  losses,  well  below  the  limit  of  detection 
of  our  apparatus,  which  is  to  say  below  about  0.002  in  loss  tangent. 

Figure  4  shows  schematically  the  principle  of  the  coupler  and  the 
electrical  arrangement  used. 

A  second  coupler  has  been  built  with  simple  cylindrical  geometry  having  a 
somewhat  longer  cavity,  but  otherwise  similar  to  the  first  one.  Both  chambers 
accommodate  cylindrical  samples  of  38-mm  diam  and  38-mm  length.  Castor  oil 
has  been  the  fluid  used  in  our  work  to  date,  and  accurate  values  for  its 
static  compressibility  are  known  [25]  over  a  wide  range  of  temperature  and 
pressure .  Some  preliminary  work  has  been  done  using  di-2  ethyl  hexyl 
sebacate,  a  plasticizer  used  by  McKinney  and  his  coworkers  at  the  National 
Bureau  of  Standards.  This  material  is,  unfortunately,  reactive  with  a  number 
of  materials  of  interest  in  our  work  and,  while  its  compressibility  is 
accurately  known,  and  while  it  has  the  advantage  over  castor  oil  of  having  a 
much  lower  viscosity  allowing  it  to  be  used  at  much  lower  temperatures,  its 
reactivity  has  ruled  out  its  further  use  in  our  work 

The  "working"  equation  for  the  coupler  has  been  given  in  Eq.  (37)  which 
relates  the  compressibilities  /Jg,  p^,  and  pQ  of  the  sample,  brass  standard, 

and  oil,  respectively,  to  the  measured  voltage  ratios  R  ,  R,  ,  and  R  .  This 
.  .  s  b  o 

equation  is  based  upon  a  calibration  of  the  system  by  substitution  in  which 
measurements  of  voltage  ratios  are  made  at  each  temperature  and  pressure  of 
interest  as  functions  of  frequency  under  three  conditions:  oil  in  the 
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coupler,  oil  plus  a  brass  sample  comparable  in  size  to  the  elastomer  sample, 
and,  finally,  oil  plus  the  elastomer  sample  itself.  Since  the  brass  sample  is 
virtually  incompressible  in  comparison  with  the  oil,  this  calibration 
procedure  actually  compares  the  compliance  of  the  oil  displaced  by  the  brass 
with  the  compliance  of  the  elastomer  sample.  A  small  correction  is  made  for 
the  compliance  of  the  brass,  but  typically  this  correction  is  of  the  order  of 
1%. 


The  choice  of  suitable  transducers  is  crucial  for  good  accuracy  and  high 
precision  in  the  measurements.  Long-term  transducer  stability  is  important 
for  maintenance  of  calibration.  It  is  also  important  that  the  transducers 
come  rapidly  to  thermal  and  electrical  equilibrium  after  changes  in 
temperature  or  pressure.  After  some  not  very  satisfactory  experience  with 
other  transducers,  we  arrived  at  a  suitable  design  using  stacks  of  samarium- 
modified  lead  titanate  disks  0.75  cm  by  0.063-cm  thick,  four  disks  to  a  stack. 
The  disks  are  made  by  the  Edo  Western  Corporation  with  silvered  faces  for 
electrical  contact.  Thin  metal  foils  are  placed  between  the  disks  to  provide 
external  electrical  leads  and  the  stack  is  cemented  together  with  epoxy.  The 
stack  is  mounted,  again  with  epoxy,  to  electrical  feedthroughs  and  surrounded 
by  an  electrically  grounded  mesh  screen  for  protection  against  accidental 
contact  with  the  sample.  The  disks  in  each  stack  are  electrically  connected 
in  parallel  giving  a  total  capacitance  of  about  1300  pF.  The  stack  resonance 
frequency  is  quite  high,  about  145  kHz.  By  comparison,  the  acoustical 
resonance  of  the  coupler  used  for  the  measurements  reported  here  falls  at 
about  8.3  kHz  with  no  sample  in  the  chamber.  Thus  the  electrical  resonance  is 
high  enough  not  to  affect  the  measurements. 

Figure  4  shows  the  early  electrical  setup  for  the  system.  In  this 
arrangement,  a  Wavetek  Model  185  sweep-frequency  generator  delivered  a 
constant  sinusoidal  voltage  of  slowly  increasing  frequency  to  the  driver 
transducer .  The  output  voltage  from  the  receiver  transducer  was  fed  into  a 
PAR  5204  lock-in  analyzer  synchronized  with  the  Wavetek  frequency.  The  PAR 
output  was  then  fed  into  a  HP  7015B  X-Y  recorder  whose  X  axis  was  driven  by  a 
signal  voltage  from  the  Wavetek  that  was  proportional  to  the  frequency  of  the 
generator.  In  this  way,  plots  were  created  of  the  receiver  transducer 
output  voltage  vs  frequency.  The  driver  input  voltage  was  held  constant  so 
the  plots  could  be  converted  into  voltage  ratios  R  for  use  in  Eq.  (37) .  The 
PAR  also  gives  quadrature  components  of  the  input  signal  that  were  recorded  in 
such  a  way  as  to  give  phase  shifts  as  functions  of  frequency.  These  phase 
shifts  were  used  to  determine  the  lossy  part  K"  of  the  bulk  modulus  whereas 
the  in-phase  voltage  ratios  determined  the  storage  modulus  K  . 
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CASTOR 

SAMPLE  OIL 


Fig.  4  -  Electrical  arrangement. 


Accuracy  and  repeatability  achieved  with  the  analog  measurement  system 
described  above  was,  at  best,  fair.  Eventually  we  were  able  to  reduce  the 
uncertainty  in  K'  to  about  ±8%  and  in  K"  to  about  ±30%.  With  this  level  of 
uncertainty  it  was  only  possible  to  determine  major  features  of  the  bulk 
modulus,  but  finer  details  and,  most  important,  the  shift  factors  could  not  be 
found  with  sufficient  accuracy  to  construct  good  master  curves.  Therefore, 
the  system  was  converted  to  a  largely  digital  measuring  setup  under  computer 
control. 

The  improved  system  uses  a  Rockland  Model  5100  frequency  synthesizer  which 
is  programmed  to  produce  a  sequence  of  fixed  frequencies  at  logarithmically 
spaced  intervals .  Control  is  provided  by  an  Apple  He  computer  working 
through  a  Rockland  Interface  Model  1488A  using  standard  IEEE  488  bus  protocol. 
The  Rockland  synthesizer  produces  a  superior  waveform  with  virtually 
unmeasurable  distortion  or  frequency  modulation.  Its  voltage  output  is  also 
much  more  stable  in  amplitude  than  that  of  the  Wavetek. 

The  receiver  transducer  output  is  still  fed  into  a  PAR  5204  lock-in 
analyzer  as  before,  but  now  all  signals  both  into  and  out  of  the  coupler  go 
through  a  HP  59307A  VHF  switch.  This  is  a  very  low-noise,  solid-state  switch 
under  computer  control  which  allows  measurements  to  be  made  of  input  and 
output  voltages  in  programmed  sequence.  The  switch  also  controls  the  outputs 
of  the  Keithley  195A  digital  voltmeters  that  read  the  PAR  in-phase  and 
quadrature  outputs.  The  DVM’s  are  under  computer  control  as  well,  as  they 
read  the  PAR  outputs  for  both  the  driver  and  receiver  transducer  voltages  and 
save  the  digitized  voltages  to  a  disk  data  file. 

Several  additional, improvements  were  made.  Better  electrical  shielding 
was  employed,  notably  in  regard  to  grounding.  Acoustical  noise  isolation  was 
provided  between  the  coupler  and  the  constant  temperature  bath  that 
thermostats  the  coupler.  The  use  of  a  much  better  bath  permitted  temperature 
control  to  a  fraction  of  a  degree.  Finally,  all  computations  are  now  made 
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from  the  raw  data  and  calibration  files  stored  on  disk,  reducing  chances  of 
computational  error  and  providing  for  flexibility  in  whatever  type  of 
smoothing  or  data  manipulation  might  be  desired.  The  data  files  are 
interfaced  to  the  NRL-USRD  VAX  computer  which  provides  parametric  curve 
fitting  routines  to  fit  experimental  K'  and  K"  values  to  the  theoretical 
master  curves.  The  net  reduction  in  uncertainty  in  the  K'  and  K"  results  has 
been  impressive.  It  is  now  routine  to  obtain  uncertainties  in  K  limited  by 
uncertainties  in  the  reported  compressibility  of  castor  oil  (about  ±2%) . 
Otherwise,  the  uncertainties  in  K'  would  be  under  *1%.  As  for  K" ,  the 
uncertainties  in  direct  measurement  are  much  better  than  before,  now  having 
been  reduced  to  about  ±8%.  However,  by  determining  K"  from  K' ,  as  discussed 
earlier,  it  is  possible  to  obtain  K"  to  approximately  the  same  accuracy  as  K1 , 
namely  about  ±1%.  At  these  levels  of  uncertainty  the  WLF  shifts  can  be 
determined  with  good  accuracy  for  most  samples  and,  from  these,  the  master 
curves  can  be  generated.  The  latter  provide  the  full  range  of  information 
about  the  bulk  modulus  and  are  the  goal  of  most  elastic  modulus  measurement 
programs . 

RESULTS 


Results  of  our  measurements  to  date  fall  into  two  categories:  (1)  low- 
precision  results  on  a  relative  wide  variety  of  rubbery  elastomers,  and  (2) 
high-precision  results  on  several  elastomer  samples  from  which  master  curves 
and  relaxation  spectra  have  been  obtained. 

The  "low-precision"  results  were  obtained  with  the  early  analog,  manually 
controlled  version  of  the  equipment.  The  uncertainty  in  K"  values  was  too 
large  to  permit  more  than  rough  qualitative  conclusions  about  the  relative 
magnitude  of  K'  and  K"  where  it  was  found  that  the  loss  tangent  (tan5=K" /K' ) 
was  small  for  all  the  elastomers  tested.  Small  in  this  case  means  less  than 
about  0.05  over  the  range  of  variables  employed,  and  usually  tan(5  was 
appreciably  less  than  0.05.  The  later  high-precision  measurements  confirmed 
these  findings  quantitatively  and,  in  fact,  showed  that  the  neoprene  samples 
had  loss  tangents  around  0.02  or  less. 

The  first  group  of  results  is  shown  in  abbreviated  form  in  Fig.  5.  The 
vertical  bars  of  the  figure  represent  ranges  of  measured  values  of  K'  at  10*C 
over  the  frequency  span  from  0.1  to  3  kHz.  K«  depends  linearly  upon  frequency 
in  this  range.  Among  the  13  different  kinds  of  elastomers  shown  in  this 
figure  are  3  natural  rubbers,  4  neoprenes,  2  polyurethanes,  and  1 
representative  sample  each  of  polybutadiene,  chloroprene  (HT1066) ,  butyl  and 
nitrile  rubbers.  1  ’ 


A  significant  point  to  note  in  Fig.  5  is  that  K'  for  all  samples  varies 

J  itttl®4.u^er  S  rans,e  of  temPeratures  and  frequencies  employed.  In 
fact,  K  falls  within  ±35%  of  the  mean  value  of  2.9  GPa,  a  value  typical  of 
orgamc  liquids.  One  can  see  that  the  neoprenes  have  generally  slightly 
higher  than  average  bulk  moduli  while  polyurethanes  and  natural  rubbers  are 
somewhat  softer  The  various  other  samples  fall  close  to  the  average.  Castor 

oi1’  SjLf11}  ilult  Whlvh  13  a  typical  organic  liquid  falls  at  about  2.5  GPa 
or  a  little  less  than  the  average  for  the  elastomers  tested.  Another  feature 

L  *eSeAeaSUrem?DtS  is  the  fact  that  the  loss  ta^ts  are  sSl  over  the 

good  i?  Wasmsufficr  rrintd’  the  accuracy  °f  the  ^  *as 

g  ,  was  sufficient  to  show  that  loss  tangents  were  generally  below  0.1. 
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Subsequent  high-precision  measurements  confirmed  that,  indeed,  they  were  less 
than  0.05  and,  ordinarily,  less  than  0.02,  again  rather  like  liquids. 


Fig.  5  -  Storage  modulus  K1  for  various  types  of  elastomers.  (Bars  give 
ranges  of  K'  values  at  10*C  for  frequencies  100  to  3000  Hz.) 


As  regards  pressure  dependence  of  the  bulk  modulus,  some  early 
measurements  were  made  at  pressures  as  high  as  20  MPa,  and  K'  was  observed  to 
increase  with  pressure,  as  it  should.  The  increase  was  scarcely  evident  at 
pressures  below  5  MPa,  however,  and  since  interest  in  the  bulk  moduli  of  the 
13  samples  measured  initially  lay  in  the  low-pressure  range  0  to  5  MPa,  the 
results  reported  here  were  all  at  2.5  MPa.  It  may  be  pointed  out  that  the 
general  formulation  of  the  theory  embodied  in  Eqs.  (4)  and  (5)  contains  the 
effects  of  pressure  implicitly  in  the  volume  term  as  expressed  in  Eq.  (5) . 

We  plan  to  explore  the  pressure  dependence  in  later  experimental  work  when  a 
form  of  Eq.  (7)  that  depends  explicitly  upon  pressure  as  well  as  upon 
temperature  will  be  used. 

Figure  6  shows  the  master  curves  for  a  typical  neoprene  sample  (5109) . 

The  reference  temperature  in  this  plot  is  27*C  which  is  the  value  of  T  at  a 

frequency  of  190  Hz  as  determined  from  a  parametric  fix  of  the  experimental 
data.  The  parameters  determined  by  the  fitting  process  are  given  in  the 
figure  (KB  is  plotted  from  the  parameters  obtained  from  the  K'  data) .  The 
experimental  points  on  the  K1  plot  are  those  measured  at  100  Hz  over  a  range 
of  T  from  -10*0  to  +40*C  at  a  pressure  of  2.5  GPa.  Other  points  measured  at 
ifferent  frequencies  from  50  to  1500  Hz  have  been  omitted  for  clarity,  but 
they  also  fall  quite  close  to  the  theoretical  line  shown.  Results  in  Fig.  6 
are  typical  of  the  high-precision,  digital  measurements  made  to  date. 
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We  note  here,  parenthetically,  that  we  have  found  a  more  convenient  way  to 
find  the  WLF  shift  parameters  than  the  one  commonly  used.  This  method  is 
treated  in  Appendix  C. 


Fig.  6  -  Master  curves  of  Kr 
for  Type  5109  neoprene. 
K^=2 . 37  GPa,  AK=1.1  GPa, 

from  the  K*  data  fit. 


and  K"  as  functions  of  logarithmic  frequency 
Parameters  for  these  curves  are  a=0.037, 
fm=190  Hz,  and  Tm=27#C.  K"  is  computed 


experiment  and  theory  agree  well,  it  is  convenient  to  characterize 
the  bulk  moduli  by  the  parameters  K;,  AK,  and  a.  From  these  three,  the  master 

curves  and  relaxation  time  spectra  can  be  constructed  to  give  the  full  picture 

of  K'  and  K"  over  all  frequencies  at  temperature  from  T  to  about  100  to  150*C 

above  Tg.  In  Table  1  below  provisional  values  of  these  parameters  are  given 

for  three  other  samples  that  have  been  measured  in  detail  by  the  high- 

of Values  lit  ^  lh4a  DOt  fiDa1’  ther  serve  to  indicate  the  ranges 

of  values  to  be  expected  for  these  parameters.  g 

in  Fig.Tis'givefbv  AK/2K' ^  ^  C°“puted  lo/s  tangent  K*/K'  at  x=0  implicit 
g  given  by  AK/2K  ,  as  can  be  seen  from  Eas.  C191  and  forrt  An 

closegtoV20%6  In  all  th*8  taDS®nt  ]?hus  comPuted  from  the  data  in  Table  1  is 
alwav.^ 1  A  ,aU  ^ 6  ®xPeriments  reported  here  the  loss  tangent  was 
sipport  ?o  ^he%r  mafltude  small«r  ^an  this.  This  lends  experimental 
dependent  666  n/reSeDted  earlfer  that  AK  is  not  constant  but 

WLF  shift  technique!^  ^  temperature’  ln  contrast  to  the  main  tenet  of  the 
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Table  1  -  Provisional  Values  of  Master  Curve  Parameters 
for  Four  Elastomers. 


SAMPLE 

K 

(GPa) 

AK 

a 

T  * 
ref 

5109  Neoprene 

2.23 

1.1 

0.037 

27 

5109S  Neoprene 

2.72 

0.81 

0.026 

19 

B252  Butyl 

2.37 

0.91 

0.033 

21 

PR1526  Urethane 

2.42 

0.83 

0.017 

18 
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Appendix  A 

HOLE  SIZE  PROBABILITY  DISTRIBUTION 


The  expression  given  by  Eq.  (1)  for  the  probability  p(q)  that  there  will 
be  exactly  q  packets  of  free  volume  each  of  size  vf  associated  with  a  given 

molecule  in  a  solid  with  N  molecules  may  be  derived  by  the  following  argument. 

Consider  N  boxes  (corresponding  to  N  molecules  in  a  sample)  to  be  filled 
at  random  by  N  balls  (packets  of  size  vf) .  The  chance  that  a  given  ball  will 

go  into  a  given  l)ox  is  1/N,  and  the  chance  that  it  will  not  go  into  a  given 
box  is  [1-(1/N)] .  Therefore,  the  chance  that  a  given  box  will  end  up  with  q 

balls  is  the  chance  (1/N)q  that  q  balls  will  go  into  the  box,  multiplied  by 
the  chance  [1- (1/N) ] (N“^  that  (N-q)  balls  will  not  go  into  the  box. 

Hpnrphe^K  a£e  N -'/(N-q)  !q!  permutations  of  the  N  boxes  with  q  balls  in  a  box. 

^  J°r  bhere  t0  be  exactlY  q  packets  of  free-volume 

Vf  associahed  with  a  given  molecule  is  the  product  of  the  three  independent 

probability  factors  above.  Since  N  is  very  large,  and  much  larger  than  q,  it 

follows  that  [1-(1/N)]  ^  is  approximately  [1-(1/N)]^. 

A  fundamental  definition  for  the  base  e  of  natural  logarithms  is 


e  1  =  lim  [1-1/N] 
N-*» 


N 


(Al) 


so  the  term  [1-(1/N)]^  becomes  simply  e 

The  permutation  factor  in  the  limit  of  large  N  becomes  approximately 

NL_  _  N(N-l)fN-2l . rN-q+1)  Nq 

qJ (N-q) !  q!  = 

and  p(q)  becomes  approximately  e_1/q! 

Using  Stirling’s  approximation  for  the  factorial  q!  gives 


(A2) 


P(q)  «  ~Z  exp  [-q(lnq-l)  -  1]  . 
i  2rq 


(A3) 
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Since  q  v^/v^  where  is  the  actual  amount  of  free  volume  associated  with 
a  molecule,  the  final  probability  that  a  molecule  has  free-volume  vL 
associated  with  it  is 


P(vf) 


_Il_  exp  (-av'/v.) 

i  2 1  1 


(A4) 


where  a^ln(q)  1  and  the  (-1)  term  in  the  exponential  has  been  dropped  since 
q>>l  in  the  present  application. 
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Appendix  B 

FUOSS  AND  KIRKWOOD  FORMULATION 


The  method  of  Fuoss  and  Kirkwood  can  be  shown  in  somewhat  abbreviated 
form  by  starting  from  expressions  for  K1  and  Kn  for  the  simple  case  of  a 
single  relaxation  time,  r.  The  equations  are  the  same  as  Eqs.  (14)  and  (15): 


K'  =  K'  +  AK  — 

1  +  (wr)2 


(Bl) 


K" 


AK 


wr 

1  +  (wr)2 


(B2) 


with  the  terms  K^,  and  AK  having  the  same  meanings  as  before. 

These  equations  may  be  generalized  easily  to  the  case  of  any  number  of 
relaxation  times 


K'  =  K'  +  AK 


l 


G. 
3=1  J 


(J^r£ 


1  +  (WTj)i 


(B3) 


K"  =  AK  I 


wr. 

L 


j=l  1  +  (ur/)d 


(B4) 


where  the  new  factors  G^  are  related  in  some  way,  as  yet  undetermined,  to  the 
distribution  of  relaxation  times,  r.  Since  the  G^  appears  in  both  K'  and  K", 

iol^ionrob^?n^%be  £eatedkm°re  eafily  in  K«,  as  is  the  case,  then  the 
solution  obtained  for  K"  can  be  transformed  for  use  in  K' . 


nnd  Kirkwood,  J.G.,  "Electrical  Properties  of  Solids.  VIII 
385  (194?)  S  ln  Pol!rvln>rl  Chloride-Diphenyl  System,'  J.  Am.  Chem.  Soc.  63, 

of  Polymers'' 3rd  Ed-  (joh" * 
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Note  that  K"  has  the  general  form  shown  in  Fig.  B1  below. 


Fig.  B1  -  Loss  modulus  K"  as  a  function  of  frequency  for  a  single 
relaxation  time.  The  abscissa  is  on  a  log  scale  where 
x=£n(w/«m)=&n(f/fm)  . 

The  subscript  m  denotes  the  value  at  which  K"  has  its  maximum.  The 
constraint  that  there  be  only  one  relaxation  time  also  means  that  there  will 
be  only  one  maximum  of  K",  or  one  inflection  point  of  K1 . 

It  is  desirable  at  this  point  to  normalize  (K*-Kf )  by  dividing  by  AK.  Let 
a  new,  reduced,  complex  bulk  modulus  be  defined  by  ° 


K*  -  K1  £  iwr.G. 

■  AK  ~  -  3  *  iB  CBS) 

J=1  J 


where  K*=K' +iK" . 

Equation  (B5)  thus  defines  a  reduced  storage  modulus,  J,  and  a  reduced 
loss  modulus,  H.  For  a  continuous  distribution  of  relaxation  times  with 
spectral  density  G(t),  Q(w)  should  have  the  following  general  form,  by 
extension  of  the  corresponding  form  given  by  Perepechko 


*Perepechko,  I. I.,  "An  introduction  to  polymer  physics,"  English  Translation 
(MIR  Publishers,  Moscow,  1981),  Ch  7. 
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«  M  ■  *■  -  dr  *  i  f-*3UL  dr 

°  o  1  +  (wr)  o  1  +  (wr)2 


(B6) 


It  remains  for  G(r)  to  be  identified  with  physical 
etc.  From  Eqs.  (B5)  and  (B6)  it  is  evident  that 


parameters  such  as  K1 , 

o' 


AK, 


J  = 


K»  -  K1 
_ c 

AK 


r”_G(r)  (wr)2  dr 
o  1  +  (wr)2 


(B7) 


/ 


G(r)  wr 
b  1  +  (wr)2 


(B8) 


The  relaxation  times  r  may  be  normalized  by  a  factor  rm  chosen  to  make  w  r  =1 

(assuming  that  H(w)  has  a  single  maximum  at  u=uj .  Then  new  variables  x^nd  s 
may  be  introduced. 


x  -  In (w/wm)  and  s  =  -ln(r/r  ) 


so 

(B9) 


wr  =  exp  (x-s) . 


the  variables  into  Eq.  (B8)  and  noting  that  in  the 
integrations  w  and  x  are  considered  fixed  so  ds=-dr/r, 


(BIO) 


Introduce  a  new  function  F(s)  defined  by 


F(s)  =  -rG(r)  =  -r  e  SG  (r  e"s) 
m  v  m  J 


(Bn) 
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Then 


H(x) 


...  f  F(s)  exp  (x-s~)  , 

J  1  +  exp  [2(x-s)j  ds 


J  F(s)  L(x-s)  ds 


(B12) 


where 


L(x) 


x 

e 


1  +  exp2x 


=  2  sech  x 


(B13) 


Note  that  Eq.  (B12)  expresses  H(x)  in  the  form  of  a  convolution  of  the 
functions  F  and  L. 

Taking  Fourier  transforms  of  both  sides  of  Eq.  (B12)  yields 


^U(H)  -  J  H(x)exp(iux)  dx  =  J  J  L(x-s)F(s)exp(iux)  dxds  = 
—to  —00  —00 


00  oq 

j[  L(t)exp(iut)  dt  J  F(s)exp(ius)  ds  = 


(L) 


FU(F) 


(B14) 


where  F^(H)  signifies  the  transform  of  the  function  H  with  transform  variable, 

u.  Solving  for  F  (F) ,  the  transform  of  F(s),  and  using  the  inversion  theorem, 
one  gets 


if"  i  r00  F  (H) 

=  2t  j[  *u  (F)exp(-ius)  du  =  |  p  exp(-ius)  du. 


(B15) 


This  equation  can  be  simplified  if  a  functional  form  can  be  found  for  the 
er?e  ’  j  ^*6  proper  form  is  given  by  Eq.  (B13) .  Fuoss  and  Kirkwood  have 
evaluated  the  Fourier  transform  of  L  by  contour  integration  giving 
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Fu(L)  "  F-u(L)  =  \  sech(fH] 


(B16) 


Therefore,  Eq.  (B15)  may  be  written 


F(s)  =  L  j  cosh (Is)  FJH) 


-1US  , 

e  ds 


r  -op 


(B17) 


Fu(H)  may  be  found  experimentally  by  taking  a  Fourier  transform  of  the  reduced 
loss  modulus  Kn(w)/AK. 


Alternatively,  let  K"/AK  [=  H(x)]  be  fitted  by  a  function  whose  form  we 
expect  to  be  a  suitably  spread  hyperbolic  secant  function  of  argument  (ax) 
where  a  is  the  spread  parameter .  In  logarithmic  coordinates  x  write 


H(x)  =  H(0)  sech(ax) 


(B18) 


as  the  fitting  function.  The  parameters  H(0)  and  a  are  to  be  determined  from 
a  best  fit  of  this  function  to  the  experimental  K"  data  on  a  logarithmic 
frequency  scale.  6 

In  th?  bulk  modulus  work,  however,  it  is  better  to  determine  a  by  fitting 
the  experimental  K' (x)  curve  rather  than  K»(x),  so  it  is  necessary  to  see  whit 

be  f°r  K'  corresP°nding  with  sech(ax)  for  K" . 
Ordinarily  this  could  be  done  most  easily  by  using  the  Kramers-Kronig 

relations  relating  real  and  imaginary  parts  of  a  linear  casual  transform  such 

In  the  present  case,  where  «<1,  there  are  mathematical  difficulties  which 
complicate  such  a  direct  determination  of  K'  .  We  therefore  proceed  along 
i  erent  lines.  Going  back  to  Eq.  (Bl)  or  (B3)  for  the  case  of  a  single 

(BirKioL^"1Vth°°ca^  a=l“Sing  losarith"ic  *”*“■“*  coordinates,  x,  Eq. 


K' (x)  =  ^  [1  +  tanh(x)] 


(B19) 


In  a  similar  way,  Eq.  (B3)  can  be  shown  to 
K  with  a=l,  namely  Eq.  (B18)  in  the  form 


give  the  Fuoss-Kirkwood  result  for 


K"(x)  =  Kn(0)  sech(x)  . 


(B20) 
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For  the  more  general  spread  case  this  becomes 


K*(x)  =  K"(0)  sech(ax) . 


(B21) 


We  therefore  take  for  the  corresponding  generalized  K' (x) 


the  expression 


K' (x)  -  +  y'  [1  +  tanh (ax)] .  (B22) 


It  should  be  noted  that  the  slope  of  the  K1 (x)  curve  given  by  Eq.  (B22)  at 
the  center  of  the  transition  region  where  x=0  is  simply  ocAK/2.  This  makes  it 
easy  to  find  the  value  of  a  from  the  experimental  data  fitted  to  the  form  of 
Eq.  (B22)  if  the  data  covers  the  inflection  region  of  K1 . 

The  pair  of  expressions  for  K'  and  K"  in  Eqs.  (B21)  and  (B22)  make  it 
possible  to  find  the  loss  modulus  from  the  storage  modulus  K'  which  is  more 
accurately  found  from  experimental  data  than  is  K" .  This  provides  a  powerful 
way  to  obtain  more  accurate  loss  modulus  values  and  is  an  important  result  of 
the  dispersion  relation  approach  to  elastic  moduli. 

the  relaxati°n  spectrum  G(r),  one  must  first  find  F(s)  using 
hq.  (B17).  This  requires  the  Fourier  transform  Fu(H)  which  can  now  be 

calculated  explicitly  from  the  expression  for  H  in  Eq.  (B18) .  Fuoss  and 
Kirkwood  have  done  this,  and  the  resulting  expression  for  G(t)  is 


G(r)  =  | 


„  r 
a  _m 

cos  ^ 

¥) 

|  cosh (as) 

r  r 

2  1 

cos 

[f 

1  2 

I  +  sinh  (as) 

(B23) 


where  s=-ln(r/rm)  as  usual. 
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Appendix  C 

DATA  SMOOTHING  PROCEDURES 


The  method  used  in  this  report  to  smooth  the  data  and  construct  a  master 
curve  will  be  described  in  detail  in  this  appendix  since  there  are  apparently 
some  differences  in  procedure  among  various  workers.  We  do  not  use  the  WLF 
method  as  described  by  Ferry  but  employ  what  may  be  considered  a  variation 
of  it  making  use  of  the  linear  form  of  the  zero-frequency  expression  for  K/T 
as  set  forth  in  Eq.  (8) ,  namely 


[In (K/T)]-1  =  d  +  bT.  (Cl) 

This  equation  permits  the  data  at  a  given  frequency  but  at  several 
temperatures  to  be  smoothed  with  respect  to  temperature  by  linear  regression 

of  [ln(K' /T) ]  1  vs  T.  This  procedure  is  carried  out  in  the  present  work  at 
two^ convenient  frequencies,  0.1  and  1  kHz  over  a  temperature  range  -10*  to 
+40  C  with  data  taken  at  5*C  intervals.  Since  K1  is  very  nearly  a  linear 
function  of  frequency  between  0.1  and  1  kHz  over  this  temperature  range, 
smoothing  at  the  two  frequencies  with  respect  to  temperature  provides  a  series 
of  short  segments  of  K1  vs  frequency  between  0.1  and  1  kHz,  one  for  each 
temperature.  These  are  plotted  on  a  logarithmic  frequency  scale,  one  above 
the  other,  labeled  by  temperature.  Then,  beginning  with  the  segment  having 
the  largest  slope,  each  segment  is  shifted  horizontally  until  it  visually  fits 
together  with  its  neighbors  forming  a  relatively  smooth  curve.  This  is 
similar  to  the  procedure  used  in  the  regular  WLF  method  to  fit  pieces  of  the 
K  (f)  curve  together.  The  fitting  process  is  easier  and  better  if  the 
segments  overlap  considerably,  and  Ferry  makes  such  overlap  a  condition  for 
proper  application  of  the  WLF  method.  However,  such  overlap  is  not  a  rigid 
requirement  when  the  following  method  is  used. 

Let  Eq.  (19)  be  written  in  three-parameter  form: 


K1 (x)  =  A  +  B  tanh  (ax) 


(C2) 


where  A=  (KUK^) /2 ,  and  B=AK/2. 

nf  US1+g  a  Par^metr^c  cufve  fitting  routine,  one  finds  the  best  values 

thl  ?ara“?ter?  A»  B'  a“d  “  the  s?read  parameter.  The  result  is  a  best  fit 
criterion110  10D  ^  Eq’  °2  ’  beSt  * lt  being  defcermined  by  a  least-squares 


Landel?  R-F-,  and  Ferry,  J.D.,  "The  temperature  dependence  of 
J.  ^  Shem  Socni?"  370^(1955) S  P°lymerS  ^  ^  «lass-forming  l^ids," 
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An  alternate  approach  is  to  perform  a  linear  regression  on  each 
experimental  K1 (f)  curve  between  0.1  and  1  kHz.  These  linear  segments  of  K1 
can  then  be  plotted  and  shifted  as  described  above,  followed  by  the  same 
parametric  curve  fitting.  This  procedure  does  not  work  as  well.  Apparently 
experimental  variations  in  the  slopes  determined  in  this  way  are  greater  than 
when  the  smoothing  with  respect  to  temperature  is  done  using  Eq.  (Cl) .  For 
this  reason,  smoothing  by  means  of  Eq.  (Cl)  has  been  used  to  obtain  all  the 
master  curves  reported  here. 

The  theoretical  slope  of  the  K1  curve  obtained  by  differentiating  Eq.  (C2) 
is 


dK1 /dx  =  a  7^  sech2 (ax) .  (C3) 

A  convenient  frequency  within  the  range  of  the  experimental  data  is 
chosen,  and  experimental  values  of  the  slopes  dK* /dx  at  this  frequency  are 
plotted  against  temperature.  The  resulting  curve  resembles  Eq.  (C3)  since  x 
is  a  function  of  temperature.  The  exact  correspondence  can  be  found  by  noting 
from  the  definitions  of  x  [Eq.  (16)]  and  the  WLF  shift  factor,  a^,,  [Eq.  (30)], 

that  if  the  shift  factor  is  referred  to  temperature  T^,  at  which  the  slope  is 

a  maximum  or  at  which  x=x  =0  on  the  log  frequency  scale,  we  have  the  simple 
relation 


x  =  -ln(a,j,) 


c  (T-T  ) 

1  v  nr 

c  +  (T-T  )  • 
2  v  nr 


(04) 


This  equation  makes  it  possible  to  translate  Eq. 
temperature  as  the  independent  variable,  giving 


(C3)  from  frequency  to 


dK1  /dx  =  a  sech2 


[c  (T-T  ) 

1  v  nr 

cn  +  (T-T  ) 

2  v  nr 


(05) 


where  c^  and  c 2  are  the  WLF  constants  referred  to  temperature  T^. 

Once  the  experimental  slopes  have  been  fitted  to  Eq.  (C3)  for  the  chosen 
frequency,  smoothed  values  of  the  slopes  can  be  used  to  improve  the  set  of 
segments  of  K1 (x) ,  one  for  each  temperature,  which  may  then  be  used  in  the 
usual  way  to  carry  out  the  WLF  shift  process  to  obtain  a  master  curve.  In  so 
doing,  a  reference  temperature  for  the  master  curve  must  be  chosen.  It  need 
not  be  T  ,  and  the  choice  is  arbitrary,  but  it  should  lie  somewhere  within  the 

experimental  temperature  range  for  best  accuracy,  as  Ferry  has  pointed  out. 
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The  experimental  master  curve  constructed  as  described  above  by  the  WLF 
method  may  then  be  fitted  to  Eqs.  (19)  or  (24)  to  find  the  parameters  K' ,  AK, 

and  a  which  fully  determine  the  shape  of  the  master  curve.  Its  position  on 
the  frequency  scale;  i.e.,  the  value  of  f  ,  is  not  determined  by  this  fit. 
however .  m 
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